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ABSTRACT 

We have investigated the radial g — r color gradients of early-type galaxies in the Sloan Digital 
Sky Survey (SDSS) DR6 in the redshift range 0.00 < z < 0.06. The majority of massive early-type 
galaxies show a negative color gradient (red-cored) as generally expected for early- type galaxies. On 
the other hand, roughly 30 per cent of the galaxies in this sample show a positive color gradient (blue- 
cored). These "blue-cored" galaxies often show strong H/3 absorption line strengths and/or emission 
line ratios that are indicative of the presence of young stellar populations. Combining the optical data 
with Galaxy Evolution Explorer (GALEX) ultraviolet photometry, we find that all blue-cored galaxies 
show UV— optical colors that can only be explained by young stellar populations. This implies that 
most of the residual star formation in early-type galaxies is centrally concentrated. Blue-cored galaxies 
are predominantly low velocity dispersion systems, and tend to live in lower density regions. A simple 
model shows that the observed positive color gradients (blue-cored) are visible only for a billion years 
after a star formation episode for the typical strength of recent star formation. The observed effective 
radius decreases and the mean surface brightness increases due to this centrally-concentrated star 
formation episode. As a result, the majority of blue-cored galaxies may lie on different regions in the 
Fundamental Plane from red-cored ellipticals. However, the position of the blue-cored galaxies on 
the Fundamental Plane cannot be solely attributed to recent star formation but require substantially 
lower velocity dispersion. Our results based on the optical data are consistent with the residual star 
formation interpretation of Yi and collaborators which was based on GALEX UV data. We conclude 
that a low-level of residual star formation persists at the centers of most of low-mass early-type 
galaxies, whereas massive ones are mostly quiescent systems with metallicity-driven red cores. 
Subject headings: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: photometry - 
galaxies: structure 



1. INTRODUCTION 

Deciphering the star formation histories of early- 
type galaxies is a key issue to understanding the fun- 
damental processes of galaxy formation and evolu- 
tion. Spatially resolved analyses can shed light on 
the build up of the stellar populations. Observations 
have shown that the colors of ea rly-type galaxies get 
bluer from the cente r outw ar ds dBoroson et al.l 11983 : 
Franx and Illingworthl Il990t IPeletier et all Il990allb : 



Michar d 1999; Idi art et all l2002t Ide Propris etal.l 12005 
Wu et alJ [20051 : lLa Barbera k. De Ca rvalho 20091. Fur 



thermore, the lack of a strong trend in the gradient over 
a wide range of lookback time ([Ferreras et al.ll2009h sug- 
gests that this negative color gradient originates from 
metallicity variations, with the center being more metal 
rich. This is often taken as evidence pointing to a sim- 
ple evolutionary model for early- type galaxies. This 
classical collapse m o d el (lEggen. Lvnden-Bell fc Sandagel 
[T^[LaJso"nlll974allbl . ll975l : lCarlberg|ll984l) suggests that 

t Einstein Fellow 



early-type galaxies form in highly efficient starbursts at 
high redshifts and have evolved without any substantial 
amount of subsequent star formation. 

While early-type galaxies were traditionally considered 
to be dynamically simple stellar systems with homo- 
geneous stellar populations, it is now clear that they 
are likely to have undergone complex and varied for- 
mation histories. Recent studies on the color gradi- 
ents of early-type galaxies have shown that a signifi- 
cant f raction feature positive color gradients (i.e. blue 
cores, Michard 1999; Menanteau, Jimenez & Matteucci 
120011 : iFerreras et all 120051 : lElmegreen et al.1 I2005D . To 
explain these positive color gradients, it is necessary 
to have age gradients in addition to metal licity gra- 
dients ([Silva fc Elstonl Il994t iMichardl l2005| ). These 
age gradients could be a natural result of the hi- 
erarchical merger scenario, where early-type galax- 
ies often form as a re sult of galaxy mergers (e.g. 
iToomre fc Toomrel 119721 ). In this model, early- type 
galaxies form as the result of successive mergers and 
are thought to have continued or episodic star for- 
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mation events. Indeed, iMenanteau. Abraham fc Ellisl 
(|2001[ ) found that a remarkably large fraction of early- 
type galaxies at high redshift exhibit positive gradi- 
ents as a result of recent m erging or an inflow of ma- 
terial. iFerreras et al.1 (|2005l ) also reported that about 
one-third of field early- type galaxies at z ~ 0.7 have 
blue cores. These findings suggest that recent episodes 
of star formation ta ke place at the center of these 
galax i es (see also IMenanteau. Jimenez fc Matteuccil 
20011: iFriaca fc Terlevichl 120011: iMarcum et all l2004t 
Elmegreen et alj|2005l) ~ 

These findings are also related to the positive gradi- 
ents found_Jnj}Oj5[^^ 

ies: IDressler fc Gunnlll983l 119921: iPoggianti et all 1 19991: 
Norton et al.ll2001l: [Bartholomew et al.ll2001l lGotoll20o4 
Yamauchi fc Goto! 120051 ). lYang et all (|2008T i suggests 
that E+A galaxies are remnants of galaxy-galaxy in- 
teractions / mergers and ultimately evolve into early-type 
galaxies (see also iGunn fc Gottl H97l IZabludoff et ad] 
H996tlBlake et al.ll2004[ ). They demonstrated that a large 
fraction of E+A galaxies have positive color gradients 
and these gradients can evolve into negative gradients as 
the populations age. 

The Galaxy Evolution Explorer (GALEX) ultraviolet 
(UV) filters allow us to detect very small amounts of 
young stars within an old stellar population. Therefore, 
GALEX data provide a unique opportunity to study the 
recent star formation history of galaxies. Many studies 
have shown that a significant fraction of early-type galax- 
ies exhibit enhanced UV light as a sign of recent star 
formation (lYi et al. 2005; S alim et al.ll2007t; [Donas et al.l 
[2007tlSchawinski et al.ll2007aHKavirai et~aLll2007l I2008E 
The observed positive color gradients in the center of 
some early-type galaxies could be naturally accounted 
for by these young stellar populations. 

The role of environment in the formation and evolution 
of galaxies is another key issue. The physical properties 
of early- t ype g alaxies correlate with their environment. 
IDressler! (|1980f ) pointed out that the environment has 
an effect on galaxy morphology, whereby the abundance 
of early- type galaxies increases in dense environments. 
This morphology-density relation implies that early-type 
galaxi es are more common in clusters. Bamford et al. 
(2009) recently suggested that the color-density relation 
is stronger than the morphology-density relation at fixed 
stellar mass. Most low stellar mass galaxies of any mor- 
phology are blue in low-density environments but red in 
high-density regions. They also found that there is a sub- 
stantial fraction of galaxies in low-density environments 
with blue colors but with an early-type morphology. It is 
also worth considering the eff ect of the nea r est ne ighbor 
galaxies on galaxy properties. iPark fc Choil (|2009t ) found 
that late-type neighbors enhance the star formation ac- 
tivity of galaxies and these effects occur within the virial 
radius. Furthermore, on a sample of close pairs only com- 
prising early- type galaxies. iRogers et al.l ([2009T ) found an 
enhancement of residual star formation with respect to 
a field sample, and a correlation between (mild) active 
galactic nuclei (AGN) activity and pair separation. 

The Fundamental Plane (FP) is a key scaling relation 
for early-type galaxies. It is a two-dimensional plane 
in the three-dimensional manifold spanned by global 
structural parameters (effective radius, surface bright- 
ness and stellar velocity dispersion; e.g. Dressier et al. 
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Fig. 1. — Our sample of SDSS DR6 visually classified early- 
type galaxies at 0.00 < z < 0.06. Green points in the background 
indicate all early-type galaxies. Orange points represent the size- 
limited sample (a criterion imposed to obtain robust color gradi- 
ents). We plot galaxies matched with GALEX as red points. We 
also plot histograms for M r and z. 

ll987HDiorgovski fc DavTslll987[ ). Through chan ges in the 
surface brightness distribution, recent star for mation can 
chang e the position of a galaxy on the FP. IChoi et all 
(2009) found that the FP of E+A galaxies is different 
from that of quiescent early- type galaxies, and this is 
most likely due to a rec ent starburst in the central re- 
gions. Ueong et al.l ([2009T ) also suggested that a dominant 
fraction of the tilt and scatter of the UV FP is due to 
the presence of young stars in galaxies with blue colors 
in the NUV, that are generally low-mass systems. 

In this paper, we measure the radial color gradients 
of early-type galaxies at 0.00 < z < 0.06 fro m Data Re- 
lease 6 of the Sloan Digital Sk y Survey (SDSS.lYork et al. 
| 2000t IStoughton et alj 120021: lAdelman-McCarthv et~ 
2008) and use them to investigate the star formation his- 
tory in early- type galaxies. Hereafter, galaxies with a 
negative color gradient are labelled as "red-cored galax- 
ies" . Similarly, galaxies with a positive gradient are la- 
belled as "blue-cored galaxies" . In § [21 we describe the 
sample and data reduction. In § [31 we present the prop- 
erties of early-type galaxies such as the color-magnitude 
relation, UV luminosities, velocity dispersion, H/3 line 
strengths, emission line diagnostics, FP, concentration 
index and environment with respect to radial color gra- 
dients. In § [?1 we quantify and discuss the recent star 
formation in early-type galaxies and revisit the effects 
of star formation on the FP. Finally, we summarize our 
findings in § [5] 

Throughout this paper we assume a ACDM cosmology 
with fi m = 0.3 and H = 70 km s _1 Mpc" 1 . 

2. SAMPLE SELECTION AND DATA REDUCTION 

2.1. Early-type Galaxy Selection in SDSS 

The Sloan Digital Sky Survey (SDSS) targets a large 
portion of the northern sky providing photometry in 
u, g, r, i and z bands. Furthermore, the SDSS spec- 
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Fig. 2. — Slope errors of early-tye galaxies versus g — r radial slope - defined as d(g — r)/dlog(R/R c ) — in the middle. The black solid 
lines indicate 0.5<r (two sloped solid lines close to the vertical dotted line) and ltr confidence level lines, respectively. We select blue-cored 
galaxies that have steeper slope than the 0.5<r confidence level line (blue stars). Red circles represent galaxies with negative g — r gradients 
(red-cored) that have steeper slope than the 0.5<r confidence level. Gradients from 0.5cr and la level confidence are plotted as open and 
filled symbols, respectively. Bigger symbols represent steeper slopes. Sample plots of radial g — r color gradients with different slope and 
slope errors are shown around the middle panel. 



troscopic survey includes virtually all extended objects 
with Petrosian ma gnitude r < 17.77 t hat escape fiber 
collision problems (jStrauss et al.l I2002D . We start by 
selecting all unsaturated galaxies with spectra in the 
redshift interval 0.00 < z < 0.06 and apply the ap- 
parent r band Petrosian magnitude cut of 17.50 to en- 
sure the imaging quality suffices for the morphological 
classification of the sampl e. The SPSS photo metry is 
given in the AB system (|Oke fc Gunnl Il983| ). How- 
ever, it is known that the photometric zero points are 
slightly offset with respect to the AB standard. There- 
fore, we use the SDSS photometric flux calibration algo- 
rithm to obtain a calibrated magnitude and apply the 
offset between the SDSS system and the AB system 
(Bohl in. Dickinson fc C alzctti 20 0~iT). We also corr ect for 
Galactic extinction using the lSchlegel et al.l (|1998l ) maps 
provided by th e SDSS pipeline and appl y a k-correction 
as described in lBlanton fc Roweisl (|2007l) . For galaxy lu- 
minosities, we used the SDSS petroMag, which is a good 
estimate of the total light, while the colors were derived 
from SDSS modelMags, which give a more reliable mea- 
sure of unbiased colors regardless of any color gradients 



(|Strauss et al.ll2002t) . 

Our next step involves the visual classification of the 
sample. We opted for a morphological classification in- 
stead of those based on color or concentration, as it re- 
duces the bias that those methods have with respect to 
the presence of young stars. We make our initial se- 
lection based on the f racDev parameter, which is the 
weight of the deVaucouleur's profile in the best compos- 
ite (deVaucouleur's + exponential) fit to the image in 
each band. We extract all galaxies whi ch have f r acDev 
> 0.95 i n g, r and i bands, following lYi et al . (12005T) 
(see also iKavirai et~afl l2007t ISchawinski et al.l booTaj ). 
On this sample, we perform visual inspection to exclude 
galaxies with spiral arms or distorted features. In ad- 
dition, we exclude small galaxies to guarantee a reliable 
measurement of the radial color gradients (see §2.2). The 
final sample comprises a total of 5,002 early- type galax- 
ies. Table [T] presents a summary of the selection criteria. 

We cross-match the detections in the GALEX Medium 
Imaging Survey (MIS) with the SDSS early-type sample 
from 704 fields, with a 4 arcsec tolerance. In the case of 
multiple galaxy matches within a given matching radius, 
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TABLE 1 
Summary of sampling criteria 



Criterion 



r < 17.50 

0.00 < z < 0.06 



Re - Rpsf > 1-188" 



Reason 
Ensure photometric depth 
Limit redshift for morphological 
classification 

Robust morphological classification 
Guarantee a reliable measurement of 
radial color gradients (having more 
than 3 fitting points between 
the PSF width and effective radius) 



we select the UV source with the smallest angular sepa- 
ration from the position of the SDSS galaxy. We also de- 
redden the colors with respect t o Galactic extinctio n us- 
ing A NUV = 8.741 x EC B-V1 (IWvder et al.l feOOSD with 
the reddening maps from lSchlegel et all (|1998l ) and apply 
a k-correction. In Figure Q] we plot all visually-inspected 
early-type galaxies as green points. The orange points 
correspond to the final size-limited sample (see §2.2 for 
details). Those galaxies in the final sample with a UV 
detection from GALEX are shown as red points. We also 
plot histograms for M r and z. 

2.2. Measurement of color gradients 

The SDSS provides the corrected frames (flat-fielded, 
sky-subtracted, and calibrated sub-images corrected for 
bad columns and cosmic rays) in five bands with 2046 x 
1489 pixels, sampled with 0."396 x 0."396 pixels. Al- 
though the images are delivered pre-processed, we per- 
form our own sky subtraction by measuring the sky level 
using the IDL "sky" Routine which adopts the tech- 
nique from the DAOPHOT because the SDSS photomet- 
ric pipeline tends to overes timate the sky background 
von der Linden et al.ll2007t ). We use the modal values 
for sky estimation. 

We carry out the surface photometry of our sample 
of early-type galaxies in g and r bands (both of which 
are redward of the Balmer break) by measuring the sur- 
face brightness along elliptical annuli, implemented in 
the ELLIPSE task within the STSDAS IS0PHQTE package 
in IRAF 2 (Image Reduction and Analysis Facility). We 
fixed the center of the isophotes to the center of the light 
distribution, and left the position angle (PA), ellipticity 
(e) and surface brightness (//) free as a function of ra- 
dius. Ellipses were fit to the higher S/N r band image 
only, and then overlaid on the g band images so that the 
colors are extracted from the same regions. 

From the radial surface brightness profiles, we derive 
g — r color gradients with a least-squares fit. The radial 
range of the fit is adjusted individually for each galaxy. 
The minimum radial coordinate of the fit is given by the 
FWHM of the point spread function (PSF) , in order to 
minimize seeing effects. The PSF widths in the g and 
r bands are different. We select the worse case between 
them. Moreover, we use in the fit pixels out to an effec- 
tive radius (Re), because the data usually become noisy 
at larger radii. From this sample, we select those galaxies 
with more than 3 points (i.e. ~1.188 arcsec) within the 
fitting range considered, to guarantee a reliable measure- 
ment of the radial color gradient - defined by its slope, = 
d(.9 — 0/dlog(R/R c ). Additionally, we perform a Monte 
Carlo simulation comprising 10 3 realizations adding noise 
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Fig. 3. — The distribution of g — r color gradients. From the 
total sample of 5,002 early-type galaxies, we select a subsample of 
761 galaxies with significant color gradients. The gray histogram 
represents the entire sample of 5,002 galaxies (most of them show- 
ing a very small color gradient). The blue solid and red dotted 
histograms represent blue- and red-cored galaxies, respectively. 

to each point compatible with the observations, to assess 
the accuracy of the color gradients measured. We plot 
the uncertainty of the slope versus color gradient in the 
middle of Figure HJ The size of the symbols is propor- 
tional to the slope. We adopt a 0.5a confidence level 
(two sloped solid lines close to the vertical dotted line) to 
pick out galaxies that show statistically-significant color 
gradients. Using this criterion, we identify hereafter the 
blue-cored galaxies in Figure [2] as blue stars. We label 
galaxies with slope < —0.5 x (slope error) as red-cored 
galaxies and show them as red circles thourghout this 
paper. Furthermore, open and filled symbols indicate 
the galaxies with slopes steeper than 0.5<r and la con- 
fidence levels, respectively. Table [2] presents the mean 
values of the slope for different subsamples as shown in 
Figure [2] Note that galaxies with smaller uncertainties 
usually have negative slopes (red-cored). Furthermore, 
for a significant fraction of the sample the error estimates 
are rather large, a result of the shallow imaging of the 
SDSS survey. 

We plot the distribution of g — r color gradients in 
Figure [3] The gray line shows the histogram for a 
total of 5,002 galaxies. The blue solid and red dot- 
ted histograms correspond to blue- and red-cored galax- 
ies, respectively. The overall fraction of blue- and red- 
cored galaxies within our criterion is roughly 4 per cent 
(216/5,002) and 11 per cent (545/5,002), respectively. 
We note that roughly 30 per cent of the galaxies with 
detectable color gradients have positive values (i.e. blue- 
cored). In other words, more early- type galaxies tend to 
have negative color gradients (i.e. red-cored). 

To measure the line strengths, we use the GANDALF 3 
(Gas AND Absorption Line Fitting) package (Sar zi et al.l 
2006). It is important to remove the contribution from 
emission lines in order to track the true values of absorp- 
tion from the underlying stellar populations. Contami- 
nation due to emission lines from extended regions of ion- 
ized gas exists in roughly 75 per cent of early-type galax- 
ies. The GANDALF measurements, such as standard 
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TABLE 2 

g — r COLOR GRADIENTS 



Sample 


mean(slope a ) 


median(slope) 


mean(Aslopc) 


Total (5,002) 


-0.009 


-0.017 


0.21 


Aslope 6 < 0.1 (294) 


-0.035 


-0.036 


0.08 


Aslope < 0.2 (1,555) 


-0.026 


-0.027 


0.13 


Aslope < 0.3 (2,564) 


-0.019 


-0.023 


0.18 


Bluc-cored c (216) 


0.25 


0.18 


0.24 


Strong blue-cored d (66) 


0.41 


0.41 


0.25 


Red-cored e (545) 


-0.14 


-0.12 


0.14 



S1 °P C - dlog(R/I 



b Slope errors of our galaxies as shown in Figure |2"1 
c Blue-cored galaxies shown in Bl, B2, B3, sBl, sB2, sB3 in Figure[2] 
d Blue-cored galaxies shown in sBl, sB2, sB3 in Figure |2"1 
e Red-cored galaxies shown in Rl, R2, R3, sRl, sR2, sR3 in Figure [2] 



Lick absorption line indices, emission line strengths and 
velocity dispersions of the SDSS galaxies are described 
in detail in Oh et al. (in preparation). 

3. RESULTS 

3.1. Color- magnitude relation 

The color-magnitude relation is widely used to 
study the star formation history of early-type galax- 
ies (s e e e.g. iBower. Lucev fc Ellisl I1992U IStanford et al.l 
119981: iFerreras et al.l I1999D . The correlation be- 
tween metallicity a nd luminosity is the ma in rea- 
son of this relation (jKodama fc Arimotol 119971 ) . The 
optical color-magnitude relation of early-type galax- 
ies reveals a sm all scatter around the mean rela- 
tion (e.g. IBower. 



mean 

,ucc\ ' fc Ellisl [1991: lEllis et a.1.1 fl997t 



3V & 

Ivan Dokkum et all 20001 ). Figure 2{a) shows the u 
color-magnitude relation. We note that the u band is 
particularly sensitive to the presence of young stellar 
populations. The figure also shows the least-squares fit 
to the whole sample (solid) and only to the red-cored 
galaxies (dotted). The gray points, blue stars and red 
circles indicate the whole sample, blue-cored and red- 
cored galaxies, respectively. A cursory inspection of this 
diagram shows that most blue-cored galaxies are located 
in the blue cloud, while red-cored galaxies reside in the 
red sequence, in agreement with re cent studies of early - 
type galaxies at moderate redshift ([Ferreras et al.l l2009). 
It should also be noted that a non-negligible fraction of 
the scatter in the optical color magnitude relation is due 
to blue-cored galaxies. 

We also show integrated u — r color versus the slope of 
the g — r color gradient in Figure 0Jb). It reveals that 
a strong correlation is present amongst the blue-cored 
galaxies, with a correlation coefficient r = 0.76. This 
implies that whenever star formation is present in early- 
type galaxies, it is mostly concentrated in the central 
regions, at least within the effective radius. We recall 
however that some galaxies show significantly bluer u — r 
colors even though they have shallower slopes. These 
galaxies show that the g — r color is significantly bluer 
out to at least the effective radius, suggesting spatially 
extended star formation. Meanwhile, we note that some 
red-cored galaxies with steep slopes show slightly blue 
u — r colors. In this case, the presence of off-centered 
young stellar populations can explain both the bluer in- 



tegrated color and the steep color gradient. 

Th e NUV color-magnitude relation (see, e.g.,lY"i et al.l 
120051 iKavirai et all 120071 : ISchawinski et al.1 l2007al ) is a 
particularly efficient tool for tracking recent star forma- 
tion, owing to its high sensitivity to young stellar pop- 
ulations. In Figure |4|c, d), we show the NUV— r color- 
magnitude relation and NUV— r color versus slope in the 
g — r ra dial gradient. T he empirical threshold at NUV— r 
< 5.5 of lYi et all (120051) (based on NGC 4552) to indicate 
recent star formation is shown as a dashed line. We find 
that blue-cored galaxies with steep slopes (filled stars) re- 
side almost exclusively below NUV— r ~ 5.5, supporting 
our findings that most of blue-cored galaxies have under- 
gone recent star formation, while red-cored galaxies are 
relatively red in the NUV. Interestingly, all blue-cored 
galaxies show blue NUV— optical colors but the converse 
is not always true. We note that there are some red- 
cored galaxies with NUV— r < 5.5. This indicates that 
some early-type galaxies have steep negative color gradi- 
ents because of the presence of young stars in the outer 
regions. For example, NGC 2974 - classified as an early- 
type galaxy in the optical - shows blue UV— optical colors 
at large radii (jJeong et al J 12007). The UV bright outer 
ring observed in this galaxy is related to a bar, demon- 
strating that recent star formation can be found in the 
form of a ring in the outer part of an early-type galaxy. 

3.2. Velocity dispersion 

The stellar velocit y dispersion is correla ted with 
a galaxy mass (e.g. iFaber fc Jackson! 119761 ). Aper- 
ture corrections are applied to the observed velocity 
dispersion due to the fact that the fixed-diameter 
fibers sample larger portions of galaxies at larger dis- 
tance, syste matically undere st imati ng the true velocity 
dispe rsi on dBernardi et all 20021: Kauff mann et all 



20031: iTremonti et al.l l2004t iKewlev et al.l 120051: 
Kewlev &: Ellison! I2008D. The appl ied aperture cor- 



rections follow ICappellari et all (|2006l ). namely: 



,Rfi 



bcrxo.066±0.035 



Cfibc 



(1) 



where Rfiber is 1-5 arcsec and R c is the effective radius. 

We show the distribution of velocity dispersions for our 
sample of early- type galaxies in the top panel of Figure [SJ 
We find that red-cored galaxies (red dotted) are rela- 
tively massive compared to average early-type galaxies 
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Fig. 4. — Gray dots in the background represent all galaxies. We plot blue-cored and red-cored galaxies with blue stars and red circles, 
respectively. A typical error bar is shown as comparison, (a) u — r vs. M r color-magnitude relation. Black solid and dotted line indicate 
the least-square fit to the total sample and only to the red-cored galaxies, respectively, (b) u — r color vs. g — r radial slope. The black line 
indicates the least-square fit of blue-cored galaxies, (c) NUV— r color-magnitude relation, (d) NUV— r color vs. g — r slope. The symbols 
are the same as Figure [2] 



(gray) while blue-cored galaxies (blue solid) tend to have 
lower stellar velocity dispersion, once mor e in agreement 
with studies at moderate re dshift (e.g. iFerreras et al.l 
2009; Ka nnappan et al.l l2009h . In the bottom panel of 
Figure OH we also plot the velocity dispersion versus 
stellar mass. The stellar masses were derived by fit- 
ting the u, g, r, i an d z SPSS photometry using a two- 
burst scenario (e.g. IFerreras fe Silkl l2000t ; iKavirai et al.l 
20071 I Schawins ki et alJl2007bl ) with the stellar models of 
Marastoiil (|2005| ). Variable parameters include the age 
and mass- fraction of the young burst, the age of the old 
burst and the amou nt of dust extinction following the 
ICalzetti et all (|2000l ) law. We measure stellar masses by 
minimizing the % 2 statistic. It is interesting to notice 
that the transition between blue- and red-cored galaxies 
occurs at around M s ~ 0.5 — 1 x 1O 1O M0, which is the 
threshold below which the stellar populations of elliptical 
galaxies become younger (jRogers et al.ll2008l ). 

3.3. Hj3 absorption line 

The presence of strong Balmer absorption lines betray 
the presence of main sequence stars formed in the past 
~ 1 Gyr. In Figure [SJ we show in the top panel the 
distribution of the equivalent width (EW) of H/3 for our 
early-type galaxies. The distribution of H/3 EWs for red- 
cored galaxies (red dotted) is similar to that of the whole 
sample (gray). Interestingly, blue-cored galaxies (blue 
solid) have a distribution of H/3 strengths that peaks at 
a roughly 1.5 A higher EW than the general sample. This 
is more evidence towards recent star formation in blue- 
cored galaxies. In the bottom panel, we show a trend 



between the H/3 EW and the g — r color gradient. We 
note that blue-cored galaxies with steeper slopes tend 
to show enhanced H/3 line strengths, which would imply 
that young stellar populations in the central region are 
responsible for the color gradients. 

3.4. Emission line diagnostics 

Emission-line diagnostic diagrams are useful tools to 
distinguish betwee n star formation and AGN activity . 
The BPT diagram (Bal dwin. Phillips, fc Terlevich|[l98lh 
allows the classification of galaxies with (narrow) emis- 
sion lines into star-forming galaxies and AGN. AGN are 
also split into Seyfert II AGN and LINERs (Low Ion- 
ization Nuclear Emission Line Region). In Figure [71 
we plot the BPT diagram using the emission line ra- 
tios [O III]/H£ and [N II] /Ha . Only galaxies with a 
signal-to-noise (S/N) above 3 (|Kauffmann et al.1 [20031 
in all four lines are s hown in the diagram. We use 
the demarcation line by Kauffma nn et all (|2003| ) (dashed 
curve) to separate star-forming galaxies. The remain- 
ing galaxies are divided into AGN and transition re- 
gion. The latter correspond to a blend of star-forming 
and AGN galaxies, and are located between the empir- 
ical line of iKauffmann et al.l (12003H and the theo retical 
maximum starburst model of iKewlev et al.l ()2001l ) (solid 
curve). AGN are subclassified into S eyferts and LIN- 
ERs as shown by the straight solid line (ISchawinsk i et all 
l2007bft . 

In Table El we show the results of the spectral line 
classification. Most early-type galaxies are virtually 
free of emission lines (~89 per cent; that is, total(100) 
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TABLE 3 

Spectral line Classification Results 



Classification 




Total 


strong blue-cored 1 


blue-cored 


red-cored 


Number 




5002 (100%) 


66 (100%) 


216 (100%) 


545 (100%) 


Strong emission 2 




549 (10.97%) 


44 (66.67%) 


85 (39.35%) 


60 (11.01%) 




Starforming... 


185 (3.70%) 


29 (43.94%) 


47 (21.76%) 


6 (1.01%) 




Transition 


176 (3.52%) 


5 (7.58%) 


12 (5.55%) 


19 (3.48%) 




Seyfert 


80 (1.60%) 


8 (12.12%) 


15 (6.94%) 


3 (0.55%) 




LINER 


107 (2.14%) 


2 (3.03%) 


11 (5.09%) 


32 (5.87%) 


Strong H/3 absorption 3 




1020 (20.39%) 


53 (80.30%) 


127 (58.79%) 


79 (14.49%) 


Quiescent 4 




3836 (76.69%) 


7 (10.61%) 


75 (34.72%) 


421 (77.25%) 



1 Blue-core galaxies which have steeper slopes than la confidence levels (filled stars). 

2 Measurement of all four lines has S/N greater than 3. 

3 H/3 EW greater than 2 A . 

4 Neither measurement of all four lines has S/N greater than 3 nor H/3 EW greater than 2 A . 




a [km s '] 

Fig. 5. — Top: The distribution of velocity dispersion for our 
early-type galaxies. The gray histogram represents the total sam- 
ple. Blue solid and red dotted histograms represent blue-cored and 
red-cored galaxies, respectively. Bottom: Stellar mass of early- 
type galaxies versus velocity dispersion. Blue-cored galaxies tend 
to have lower stellar velocity dispersion, while red-cored galaxies 
are more massive. The symbols are the same as in Figure [2] 

minus strong emission(10.97)). Twenty three per cent 
of total early-type galaxies (that is, total(lOO) minus 
quiescent(76.69)) show significant emission lines (S/N > 
3 in all four BPT lines) or strong H/3 absorption lines 
(EW{H/3) > 2). We find that about 65 per cent (that 
is, total(100) minus quiescent(34.72)) of blue-cored 
galaxies show strong H/3 absorption and/or emission 
lines. Most of the emission line ratios are consistent with 
ongoing star formation. Only 10 per cent (quiescent; 
10.61 per cent) of blue-cored galaxies with slopes steeper 
than the la confidence level (i.e. strong positive color 
gradient galaxies) show neither emission nor strong 
Balmer absorption lines. Interestingly, some blue-cored 
galaxies are classified as Seyferts. Therefore, we note 




1 2 3 4 5 

H EW [A] 



Fig. 6. — Top: The distribution of H/3 absorption line strength 
for our early- type galaxies. The gray histogram represents the 
total sample, blue solid histogram represents blue-cored galaxies 
and red dotted histogram represents red-cored galaxies. Bottom: 
g — r radial slope of early-type galaxies versus H/3 absorption line 
strength. Blue-cored galaxies show stronger H/3 absorption than 
red-cored galaxies. The symbols are the same as in Figure [2] 

that blue-cored galaxies are mainly due to centrally 
concentrated star formation and a smaller fraction could 
be associated with Seyfert activity. Most red-cored 
galaxies, on the other hand, are classified as quiescent 
(78.98 per cent). If red cores show emission lines, most 
of their emission line ratios are consistent with LINER 
AGN. red-cored galaxies appear not to be host to active 
star formation or powerful AGN. 

3.5. Fundamental Plane 

Early-type galaxies exhibit a scaling relation be- 
tween the effective radius i? c , the effective mean sur- 
face brightness fi e and central velocity dispersion a 
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Fi g. 7. — The emission line d iagnostic diagram for our sam- 
ple (Baldwin. Phillips, & Tcrlevich 1981). We only plot galaxies 
where [N II], Ha, [O III] and H/3 lines are detected with S/N > 3. 
T he curve labeled Ka03 is the empirical purely star- forming llimit 
of IKauffmann et al.1 Q2003T I. while the curve labele d KeOl repre- 
sents the theoretical maximum starburst model from ( Kcwl ev et alj 
120011) . All galaxies between the Ka03 and the KeOl are transition 
systems that they do have emission lines with c ontributions from 
both s tar-formation and AGN. The line from S07 ( Schawins ki et all 
I2007bf ) divides Seyferts and LINERS. The symbols are the same as 
Figure [2] 
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Fig. 8. — The Fundamental Plane of early-type galaxies, blue- 
cored and red-cored galaxies are marked with blue stars and red 
circles, respectively. Gray points in the background indicate all 
early-type galaxies. There are indications that red- and blue-cored 
galaxies may lie on different planes. The symbols are the same as 
m Figure [2] 

(i.e. the Fundamental P lane, FP, iDressler et al.l 119871 : 
iDjorgovski fc Davis! Il987t ). We compare here the loca- 
tion of blue-cored and red-cored galaxies on the FP. In 
order to construct the fundamental plane of our early- 
type galaxies, we have to apply a number of corrections 
as follows. We conve rt an elliptical apert ure to an ef- 
fective circular radius ([Bernardi et al.ll2003l) and apply a 



k-correction and the cosmological dimming effect: 

R = v/b^Rdev, (2) 

fi = m dcv + 2.51og(2^(R c ) 2 ) - K(z) - 101og(z + 1), (3) 

where b/a is the ratio of the minor and major axes, 
i?dcv and mdcv are the effective radius and apparent 
magnitude, respectiv ely, from the de Vaucouleurs fit 
(|de Vaucouleurslll948l ). 

We show in Figure[S]the FP of the total sample (gray), 
blue-cored galaxies (blue stars) and red-cored galaxies 
(red circles). Red- and blue-cored galaxies appear to lie 
on different planes. A more detailed analysis is on the 
way. We will discuss this discrepancy via stellar popula- 
tion modelling in §4.2. 

3.6. Concentration index 

The concentration index is often used as a quan- 
titative measure of g,alaxy morphol ogy (see e.g. 
iMorgan fc Maya!] Il957t IConselicel 12006) and is a use- 
ful tool because the surface brightness distribution of a 
galaxy depends strongly on its formation history. The 
concentration index, C, is denned as the ratio between 
Rgo and R50 from the SDSS pipeline parameters, which 
are the radii that encompass 90 and 50 per cent, respec- 
tively, of the Petrosian flux, namely: 



where x =u, g, r labels the SDSS band used t o deter- 
mine the concentration ()Stoughton et al.l I2002T) . After 
unresolved sources (which just result in a PSF), early- 
type galaxies are the mo st concentrated syst ems with 
C r > 2.6 (see Figure 1 of iFerreras et al1[2005l ). 

In Figure [SI we show the u, g and r band concentration 
indices versus the slope of g — r color for our sample of 
early-type galaxies. The color code for the symbol is the 
same as in Figure 2. The size of the symbols is propor- 
tional to the H/3 line strength. The right panel indicates 
that blue-cored and red-cored galaxies show similar con- 
centration index in the r band. However, the other pan- 
els show that a significant fraction of blue-cored galaxies 
are more concentrated in u and g bands, while there is a 
tendency for red-cored galaxies to be less concentrated, 
although examples of galaxies showing the opposite be- 
havior can also be identified. We note that the most con- 
centrated blue-cored galaxies in u band tend to also show 
strong H/3 absorption lines, indicating large fractions of 
young populations in the center. This difference as well 
as the increase in scatter are naturally explained by the 
different sensitivity of the passbands to young stars. 

3.7. Environment 

We find evidence that a large fraction of blue-core 
galaxies have undergone recent star formation from H/3 
absorption line strengths and emission line diagnostics. 
We now address the possible cause of this recent star for- 
mation in early-type galaxies. One possible scenario in- 
volves galaxy mergers and interactions. If galaxies expe- 
rience merging events involving gas-rich late-type galax- 
ies, it might result in enhanced star formation during 
the encounter. We investigate this possibility by look- 
ing at environment and also by checking the nature and 
location of the nearest neighbor. 
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Fig. 9. — The u, g and r band concentration indices of early-type galaxies as a function of g — r slope. Gray dots in the background 
indicate the total sample of visually classified early- type galaxies. The symbols are the same as Figure [2] but the symbol size represents 
H/3 absorption line strength. The typical 1-sigma error in the measurements is shown in each panel. The slope errors are measured in this 
study, and the concentration index errors are adopted from the SDSS database. 



1.0 



0.8 



-o 



0.6 



0.4 



0.2 



0.0 




KS test (P = 0.06) 



log(p 8 ) 

Fig. 10. — Cumulative distribution of the local density log(p 9 ). 
In order to minimize a selection bias with stellar mass, this figure 
shows both blue-cored and red-cored galaxies in the same stellar 
mass bin (10.5 < log(M/MQ) < 12). The KS test discriminates be- 
tween these two distributions at a 95% confidence level, blue-cored 
galaxies tend to prefer slightly lower density regions compared to 
red-cored galaxies. 
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Fig. 11. — Fraction of pair galaxies with respect to (visually clas- 
sified) morphology. The morphologies of pair galaxies are divided 
into three types: early- type galaxy (E), dwarf galaxy (Dwarf) and 
late- type galaxy (L). Blue-cored galaxies appear more are likely to 
have a late-type companion. 
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Fig. 12. — The distribution of separation in close pairs (within 
300 kpc) comprising an early- and a late-type galaxy. Pairs with a 
blue-cored galaxy are found at shorter separations with respect to 
pairs with a red-cored galaxy. 

Figure [10] shows the cumulative di stribution of the 
local d ensity par ameter p a define d by ISchawinsk i et al.1 
(|2007aD (see also lYoon et~alll2008[ ). This density param- 
eter (p g ) represents a weighted local number density, and 
is obtained by counting all neighbors within 2 Mpc with 
a Gaussian-weighting scheme. In order to avoid selection 
bias with respect to mass, we limit both blue-cored and 
red-cored galaxies to the same stellar mass range between 
10.5 < log(M/Mo) < 12. We perform the Kolmogorov- 
Smirnov test (KS test) and find that blue-cored and red- 
cored galaxies are not consistent with being drawn from 
the same parent population at the 95 per cent confi- 
dence level. These two cumulative distributions of p g 
show that blue-cored galaxies prefer lower density envi- 
ronments compared to red-cored galaxies. However, this 
may also come from the fact that more massive galaxies 
reside in denser regions. For a narrower mass bin, we 
find less clear density dependence. But we had to use 
such as large mass range for this statistical test because 
the number of blue-cored galaxies in our sample is lim- 
ited. We would need a much larger sample to break the 
mass-density degeneracy. 

We also consider the nearest neighbor galaxies and 
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TABLE 4 
Pair galaxy fraction 

Classification Fraction (r p < 120 kpc) 
Total 764/5002 (15.27%) 
blue-cored 45/262 (17.18%) 

rcd-corcd 120/647 (18.54%) 
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Fig. 13. — The mean radial g — r color gradients of blue-cored 
(blue solid) and red-cored (red dotted) galaxies. The shading cor- 
responds to the scatter within the sample, given at a la confidence 
level. 

their morphology. Two galaxies are considered to be in 
a close pair if r p < 120 kpc (where r p means projected 
physical separation) and Az < 0.001. We show the frac- 
tion of galaxies in close pairs in Table |3J Radial color 
gradients seem insensitive to pair galaxy fraction. How- 
ever, this analysis is hampered by projection effects that 
can mimic true pairs which are more serious in denser 
regions. 

If the recent star formation in the central region of 
blue-cored galaxies is caused by gas inflows from a neigh- 
boring galaxy, one might wonder whether the morphol- 
ogy of a companion galaxy - which correlates with gas 
fraction - could affect color gradients. To test this point, 
we perform visual classification of the companion galax- 
ies. The morphology of the companion galaxies are clas- 
sified into three types: early-type (E), dwarf (Dwarf) and 
late- type (L). Figure [TT] shows the fraction of pair galax- 
ies with respect to the companion morphology. Taken at 
face value, blue-cored galaxies have more late-type com- 
panions than red-cored galaxies. 

Figure [T^] shows the distribution of pair separation in 
systems comprising an early- and a late-type galaxy. We 
find that pairs with a blue-cored galaxy have smaller sep- 
arations than those with a red-cored galaxy, in agreement 
with the picture that a nearby, gas-rich late-type galaxy 
can contribute to the onset of residual star formation. 
However, we should emphasize here that most (92 per 
cent) of the blue-cored galaxies in our sample do not 
have a late-type companion. The presence of a late-type 
companion may affect the star formation of the galaxy in 
question, but not all star-forming early types are caused 
by the presence of a late-type companion. Hence, other 
mechanisms (such as minor mergers with nearby gas rich 
dwarf galaxies or even gas clouds) should be invoked to 
explain the general trend. 

4. DISCUSSION 
4.1. Age and mass fraction of young stars 
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Fig. 14. — The normalized likelihood distribution in age and mass 
fraction for the young stellar component, (a) The star symbols 
represent the best fits to the blue-cored galaxies (as classified in 
Figure 2) along with error bars in grey. Only those with x 2 <j < 3 
are shown here. The contour gives the la confidence level from 
a Monte Carlo simulation (see text for details), (b) Normalized 
age distribution, (c) Normalized mass fraction of the young stellar 
component, (d) A typical x 2 marginalization for a galaxy. This 
particular case is for the circled galaxy in (a). 

The aim of this section is to quantify the amount of 
recent star formation in our sample. In order to study 
the radial distribution of young stars, we construct the 
mean g — r radial color gradient of the sample, shown 
in Figure [T31 The dotted and solid lines indicate red- 
cored and blue-cored galaxies, respectively. The shading 
in each case represents the la confidence level. Color gra- 
dients in early-type galaxies are usually a result of the 
underlying stellar populations, mainly in age and metal- 
licity. It is difficult, however, to distinguish the effects 
of a small change in age from those of a small change 
in metallicity. Most red-cored galaxies, however, tend 
to show overall red u — r colors (see §3.1) and weak H/3 
line strengths (see §3.3), suggesting an old population. 
Therefore, we consider the negative slopes of red-cored 
galaxies as a result of a simple metallicity gradient and 
assume that blue-cored galaxies have the same metallic- 
ity gradient. This assumption is also motivated by the 
lack of evolution in the color gradient of red -cored galax- 
ies with redshift (|Ferreras et al.ll2005l[2009t) . 

To calculate the characteristic metallicity of each re- 
gion in the galaxy, we assume a single uniform age of 
12 Gyr but change the metallicity distribution until it 
matches the observed g — r co lor. The base model in this 
study are those of [YJ (l2003[ ) which are specialized for 
old populations. Since those models do not cover ages 
younger than 1 Gyr, however, w e com bine them with 
the models of iBruzual fc Chariot! (|2003| ). According to 
the models, the central region in red-cored galaxies has 
a metallicity distribution with a peak around the solar 
value, while at the effective radius the peak decreases to 
half the solar. 

We consider a two-stage star formation history to de- 
termine the age and mass fraction of the young stellar 
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Fig. 15. — The evolution of the color gradient for blue-cored galaxies. The blue solid and red dotted lines represent the mean radial 
g — r color gradients of blue-cored and red-cored galaxies, respectively. The open circles correspond to models using the observations of 
red-cored and blue-cored galaxies as constraints on a two- burst scenario (see text for details). As the populations age (from left to right), 
the blue-cored galaxy transforms into a red-cored galaxy. 
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Fig. 16. — Evolution of slope obtained from the probability 
distribution in the age and mass fraction of the young component. 
We show how the slope evolves for various possible choices of the 
young stellar mass fraction. We choose mass fractions from the 
peak of the probability distribution (solid), or at the la (dotted), 
2<r (dashed) and 3<r (dot-dashed) confidence levels. The slope of 
blue-cored galaxies decreases rapidly from positive to negative in 
about 1.5 Gyr after the star formation episode. 

component of blue-cored galaxies. This model combines 
an old and a young stellar population (iFerreras fe Silkl 
l2000t iKavirai et all (20071: iSchawinski et al.ll2007bD . The 
old population has a fixed age of 12 Gyr with a metallic- 
ity gradient determined from the mean radial color gradi- 
ent of red-cored galaxies. The young component has the 
solar metallicity and is allowed to vary in age (0.01 < 
tyoung <10 Gyr) and mass fraction (10~ 4 < fyoung < 1)- 
To constrain these two parameters, we fit the observed 
g — r color and H/3 EW to models and compute the as- 
sociated x 2 statistic to obtain a probability distribution 
of the age and mass fraction of the young stellar pop- 
ulation. We further perform a Monte Carlo simulation 
with 10 3 realizations to explore the parameter space for 
minima and the associated errors around them. The best 
fit and confidence levels are shown in Figure HUJa) and 
also in the accompanying histograms (b and c) . In panel 
(d) we show the \ 2 minimization performed on a sam- 
ple galaxy, circled star in (a). These clearly show that 
young (~ 10 2 Myr old) stellar populations are present 
in the central region and the stellar mass fraction of the 
young component reaches up to a couple of per cent. 

By evolving the young component in time, we can de- 
rive the evolution of the slope. Figure [TBI shows the evo- 



lution of model blue-cored galaxies, assuming the mean 
properties and scatters for the young population derived 
here. The color gradients observed suggest that this 
young component should have a decreasing contribution 
with increasing radius. Hence, we choose a central value 
for the mass fraction in young stars of 0.7 per cent and 
decrease this mass fraction outward in a manner that 
reproduces the observed g — r profiles. Its mean age is 
assumed to be 100 Myr to begin with, and we monitor the 
color gradient as time goes on. The solid and dashed lines 
represent the mean radial g — r color gradients of blue- 
cored and red-cored galaxies, respectively. The color 
gradients caused by the young population revert to the 
fiducial values of red-cored galaxies after about 1 Gyr, 
as expected given the small fractions in young stars. In 
Figure [16] we show the evolution of the change in slope 
with respect to different values of the mass fraction in 
young stars. We choose four different cases based on the 
probability distribution of the modeling (see Figure I14|) . 
We consider the peak of the probability distribution (0.7 
per cent, solid), and the la (1.8 per cent, dotted), 2er 
(3.9 per cent, dashed), and 3<r (5.7 per cent, dot-dashed) 
confidence levels. The slope of the blue-cored galaxies 
changes rapidly from positive to negative. The positive 
gradients in early-type galaxies must be transient fea- 
tures that are visible only for ^0.5 - ~1.3 billion years or 
so after a centr ally-concentrated sta r formation episode 
(sec figure 8 of iFerreras et aL1l2009h . 

4.2. The effect of young stars 

We now discuss the effect of the young stellar popula- 
tions on the concentration index and the FP. To explore 
the effect of young stars on these properties, we com- 
pute the change of the surface brightness distribution of 
a galaxy due to recent star formation in the central re- 
gion. We start from the Sersic profile: 

I(R) =I(0)exp[-k(R/R e ) 1 / n ] (5) 

For n > 1, k satisfies the relation k — In — 0.324 (jCiottil 
[1991 . 

We assume that the intensity profile of the underlying 
old population has the metallicity gradient constrained 
by the color gradient and that the young stellar pop- 
ulations are more centrally concentrated to satisfy the 
observed positive color gradients. A model effective ra- 
dius of 5 arcsec is chosen for this exercise. We construct 
the total light profile as the superposition of the sur- 
face brightness profile for the old and young populations. 
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Fig. 17. — Top: Concentration indicies of u, g and r band as a function of Sersic index. The red dotted line represents the concentration 
index assuming a purely old population. The blue solid line shows the effect of a centrally concentrated young population on each passband. 
Bottom: Concentration indicies versus slope. The stars represents observed blue-cored galaxies and the squares show the results from the 
modelling (see text for details). The gray horizontal line marks the mean concentration index of the total sample in each passband. 



This profile allows us to quantify the effect of a young 
component on the effective radius, mean surface bright- 
ness or concentration index. 

4.2.1. Concentration index 

The concentration index is a useful tool because 
the surface brightness distribution of a galaxy depends 
strongly on its formation history. Star formation, merg- 
ers and accretion events might affect the surface bright- 
ness distribution by rearranging the stellar components 
in the progenitors as well as by forming new stars in ge- 
ographically biased ways. For example, if star formation 
is triggered towards the center, the concentration index 
will increase. 

Figure [17] shows the concentration index of the old 
(dotted) and the mixed (solid) populations. The top 
panels of Figure [17] show that the strongest effect is 
on the concentration measured in the u band because 
of its highest sensitivity to young stars. We also show 
the increase of the concentration index as a function of 
slope in the bottom panels. The horizontal line indi- 
cates the mean concentration index for the total sample 
in each passband. Using the estimated best fit age and 
mass fraction of the young stellar component of blue- 
cored galaxies (see §4.1), we compare our models (black 
squares) with observations (stars). The concentration 
index in the model shows a clear correlation with the ra- 
dial color gradient slope. The correlation is stronger for 
shorter wavelengths. The match seems poor for the it 
band. This is probably because of the inaccuracy of our 
model assumptions. We assumed that recent star forma- 
tion occurred only in the central region, from center to ef- 
fective radius. However, the SDSS concentration indices 
are measured from the whole radial range. If the resid- 



ual star formation extends outside the effective radius, 
it would lower our model concentration index to match 
the data better. But such a detailed modeling is beyond 
the scope of this investigation. Despite the simplicity 
of the model considered here, we successfully match the 
observed concentration indices in various bandpasses. 

4.2.2. Fundamental plane 

In order to examine the effect of the young population 
on the FP, we calculate the change of the effective radius 
and the mean surface brightness due to recent star for- 
mation in the central region. We assume that the surface 
brightness distribution of a red-cored early-type galaxy 
follows a de Vaucouleurs' profile (i.e. Sersic index n — 
4). 

In Figure IToT a). we show the intensity profile of the un- 
derlying old population (dotted) and the combined pro- 
file (solid) . We show the cumulative luminosity profiles in 
Figure fTSTb). A population with young stars at the cen- 
ter tends to have a smaller effective radius and a brighter 
mean surface brightness, as shown for the Kormendy re- 
lation in Figure ITSf c). In order to study the effect of 
this centrally concentrated young population on the FP, 
we compare our models with observations, as shown in 
Figure ITSl' d) . This suggests that although the models 
predict a significant offset in the Kormendy relation, the 
effect is more subtle on the FP. Hence, the clear differ- 
ence between the observed blue- and red-cored galaxies 
on the FP (Figured]) seems beyond what can be modeled 
by a simple addition of young stars. Instead, one needs to 
consider that blue-cored galaxies indeed have on average 
lower velocity dispersions than the general sample. 

To further illustrate this point, Figure Q3J] shows the 
effect of centrally concentrated recent star formation on 
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Fig. 18. — (a) The surface brightness distribution. The dotted line represents the intensity profile assuming a purely old population. The 
intensity profile is normalized to 1(0). The solid line shows the combined light profile of the old and young populations, (b) Cumulative 
light profile. Dotted and solid lines represent the underlying old population and the combined model, respectively. The vertical lines 
represent the effective radius of each profile, (c) The effect of a centrally concentrated young population on the Kormendy relation. Red 
circle and blue star represent old and total (old plus young) model, respectively. The gray dots in the background indicate all observed 
early-type galaxies, (d) The effect of a centrally-concentrated young population on the FP. Symbols are the same as in (c). 




1.0 1.5 2.0 2.5 1.0 1.5 2.0 2.5 

log a + 0.2(/z -20.09) log a + 0.2(/z -20.09) 



Fig. 19. — The effect of a centrally concentrated episode of star formation on the FP. Blue solid and red dotted lines delimit the observed 
FP (shown at the la confidence level) for blue-cored and red-cored galaxies, respectively. The vectors show the contribution of a 100 Myr 
young components with a mass fraction of 0.7% or 5.7%. Left: Effect of a young component on red-cored galaxies. Right: Effect of ageing 
of the young component of blue-cored galaxies. 
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the FP. The solid and dotted contour lines delimit the 
region occupied by the FP for blue-cored and red-cored 
galaxies, respectively. The black and purple/pink ar- 
rows correspond to the presence of 0.7 or 5.7 per cent 
of a 100 Myr population at the center. The left panel 
considers what would happen if this young component 
is added to the population of red cores. On the right 
panel, we consider the fading of this young component 
in blue-cored galaxies, as the galaxies age. In both cases 
we find that blue/red-cored galaxies belong to different 
distributions, separated by velocity dispersion (or stellar 
mass). One can clearly see that although the most mas- 
sive blue-cored galaxies could be explained by a recent 
episode of star formation (comprising 6 per cent of its 
stellar mass content, which is hardly unjustifiable), for 
the majority of (low mass) blue-cored early-type galax- 
ies, one must consider that they correspond to systems 
with markedly lower velocity dispersion than red-cored 
galaxies. 

5. SUMMARY 

We present radial g—r color gradients of visually classi- 
fied early-type galaxies at 0.00 < z < 0.06 selected from 
Data Release 6 of the Sloan Digital Sky Survey. Due 
to the shallow exposure of DR6 color gradients are diffi- 
cult to measure in the bulk of the sample. The brightest 
galaxies exhibit a negative color gradient (centrally red- 
der) as expected for their observed radial metallicity vari- 
ations. However, a significant fraction - roughly 30 per 
cent which might depend on data selection strategy - 
features blue cores (positive color gradients). Evidence 
is presented, suggesting that these blue-cored galaxies 
have undergone a recent episode of star formation that is 
centrally concentrated. Most are blue also in terms of in- 
tegrated u — r colors and show strong central H/3 absorp- 
tion line strengths and/or emission ratios that are con- 
sistent with star forming populations. Combining SDSS 
and GALEX UV photometry, we find that all blue-cored 
galaxies show blue UV— optical colors. They also tend 
to have a lower stellar velocity dispersion. We inves- 
tigate their environmental dependence and star forma- 
tion properties. Blue-cored galaxies exhibit a tendency 
to live in lower density environments compared to red- 
cored galaxies and, if found in close pairs, have a pref- 
erence for a late-type companion. On the other hand, 
red-cored galaxies, which are relatively massive, are lo- 
cated on the red sequence as shown in the integrated u— r 
color-magnitude relation, and those with emission lines 
reside in the LINER region of a BPT diagram. 

Based on a simple model that overlays a young stel- 
lar component over an old population, our sample gives 
a mass fraction in young stars below 2 per cent with 
age 100—800 Myr in the central regions of blue-cored 
galaxies. Furthermore, these positive color gradients in 
early-type galaxies are visible only for ~0.5 - ~1.3 bil- 
lion years after a star formation event. Although this 
centrally located star formation episode can decrease the 



effective radius and increase the mean surface brightness, 
the FP of the (majority of) low mass blue-cored galax- 
ies cannot be reproduced by any amount of recent star 
formation on a red-cored galaxy. Instead, they require a 
lower velocity dispersion, i.e. blue-cored and red-cored 
early-type galaxies are, perhaps fundamentally, different 
populations, mainly split with respect to stellar mass or 
velocity dispersion. 

From the perspective of our investigation on the op- 
tical color gradients, we conclude that star-forming 
early types are not simply the star-forming counter- 
parts of the quiescent ones. While many parameters 
are connected to each other in degenerate manners, 
the most important role seems to be played by galaxy 
mass. This star-formation dichotomy (between the star- 
fo rming and quiesce n t early types) is consistent with that 
of lSchawinski et all ((2006) and also with the kinematic 
dichoto my pr e viously addressed (e.g . Davies et alj|1983t 
Bender! [l98llKorinendv fc Bender) Il996t iFaber et al 



19971 iRest et aljfeOQll lLauer et all [20071 iKimm fc Yi 
2007). Recent theoretical models fe.g. iKhochfar fe Silk 
20091 ) explain the star-forming dichotomy in connection 



with the kinematic dichotomy through distinct merg- 
ing histories by mass. More massive early types are 
interpreted as results of equal-mass dry mergers. This 
is compatible to t he axis ratio distribution analysis of 
IKimm fc Yil (|2007l ). On the other hand, small early types 
are suggested to be a product of early-type and late-type 
merger, in which case the cold gas from the late-type 
component would linger in the me rger remnant and get 
used for residual star formation. iKavirai et al.l (2009) 
further point out that it was more likely minor merg- 
ers than major ones that was involved in the formation 
of low-mass early types exhibiting signs of residual star 
formation. All these works together seem to suggest a 
coherent picture for the formation of early-type galaxies. 
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